THE MAXIMAL SPECTRAL RADIUS OF A DIGRAPH 
WITH (m + 1)2 - s EDGES 



JAN SNELLMAN 

Abstract. It is known that the spectral radius of a digraph with 
k edges is < Vk, and that this inequahty is strict except when k is 
a perfect square. For k = m^+i, £ fixed, m large, Friedland showed 
that the optimal digraph is obtained from the complete digraph on 
m vertices by adding one extra vertex, a corresponding loop, and 
then connecting it to the first [^/2J vertices by pairs of directed 
edges (for even i we add one extra edge to the new vertex). 

Using a combinatorial reciprocity theorem by Gessel, and a clas- 
sification by Backelin on the digraphs on s edges having a maximal 
number of walks of length two, we obtain the following result: for 
fixed < s 7^ 4, = (m + 1)^ — s, m large, the maximal spectral 
radius of a digraph with k edges is obtained by the digraph which 
is constructed from the complete digraph on m + 1 vertices by 
removing the loop at the last vertex together with [s/2j pairs of 
directed edges that connect to the last vertex (if s is even, remove 
an extra edge connecting to the last vertex). 

Resume. On salt que le rayon spectral d'un graphe oriente 
avec k arcs est < a/A:, et que cette inegalite est stricte sauf quand 
k est un carre parfait. Pour k = m? + i, £ fixe, m grand. Fried- 
land a prouve que le graphe oriente optimal est obtenu a partir 
du graphe oriente complet a m sommets en ajoutant un sommet 
supplementaire, une boucle correspondante, et en le reliant aux 
premiers [£/2\ sommets par des paires d'arcs. 

En utilisant un theoreme combinatoire de reciprocite de Ges- 
sel, et une classification due a Backelin des graphes orientes a 
s arcs ayant un nombre maximal de chemins de longueur deux, 
nous obtenons le resultat suivant: pour s ^ 4 fixe, < s, 
k = (m -|- 1)^ — s, m grand, le rayon spectral maximal d'un graphe 
oriente a k arcs est obtenu pour le graphe oriente construit a partir 
du graphe oriente complet a m -I- 1 sommets en enlevant la boucle 
d'un sommet quelconque ainsi que les [s/2j paires d'arcs reliee a 
ce sommet. 
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1. Introduction 



By a digraph we understand a finite directed graph with no mul- 
tiple edges, but possibly loops. Let G{m,p,q) be the digraph on 
{1, . . . , m + 1}, where there is an edge from i to j ii i, j < m or ifi < p 
and j = m+1 or if 2 = m+1 and j < q. Let M{G{m, p,q)) = M{m, p, q) 
denote the adjacency matrix of G{m,p,q); it is a 0-1 matrix with 



m' 



+ P + q ones. If Ia,b denotes the ax b matrix with all ones, and 



the matrix with all zeroes, then 



M{m,p, q) 







Il.ij Ol,m-9 






(1) 



For < £ < 2m + 1, we put M{m,£) = M{m, \i/2]J - \i/2]), i.e. 
the M{m,p, g)-matrix with p + q = i,p > q, p — q minimal. We denote 
the corresponding digraph with G{m,i). As an example, G(5,2, 1) = 
G{5, 3) is shown in figure ^ and has adjacency matrix 

■1 1 1 1 1 r 



M(5,2,l) = M(5,3) 



(2) 




Figure 1. G(5,2, 1) = ^(5,3) 

Friedland |Fri85j showed that for a fixed i there is an L{i) such that 
for m > L{i) the maximal spectral radius of a 0-1 matrix with -|- i 
ones is achieved by the matrix M{m,i). He conjectured that for any 
m,i, < i < 2m + 1, the maximal spectral radius of a 0-1 matrix 
with -|- i ones is achieved by some {m + 1) x [m + 1) matrix. For 
+ i = [m + 1)'^ — 4, he showed that the optimal matrix is not Mm,e 
but rather 



Im— l.m— 1 


Im-1,2 


^2,m-l 


©2,2 



(3) 
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It is reasonable to believe that this is the only exception and that for 
other i, M{m, i) is optimal. 

We show a weaker result: for a fixed s 7^ 4 there is an S{s) such that 
for m > S{s) the maximal spectral radius of a digraph with (m + 1)^ — s 
edges is achieved by the digraph G{m, 2m + 1 — s). 

Our main tools are the following: 

(1) A combinatorial reciprocity theorem by Gessel |(7es77j which 
asserts that for a digraph A, the generating series 

00 

HA{t) = Y,Xn{Ar, (4) 

n=0 

of Xn{A), the number of walks of length n — 1 in A, is related to 
the series of the complementary digraph by t) = 1, 

(2) A classification by Backelin |Bac95j of the digraphs of s edges 
with maximal number of walks of length 2. 

The proof runs as follows: Backelin's classification shows that for s > 6, 
m sufficiently large, the digraph with the following adjacency matrix 
has the maximal number of walks of length 2 among digraphs with s 
edges and m + 1 vertices: 





'Dlm-s/2,1 




Is/2.1 


Il,m-s/2 


1 



(5) 



The generating series for walks in that graph is 

1 + (m + l)t + + ct^ + O(t^), 

so the generating series for the complementary graph, which has adja- 
cency matrix M{m, 2m + 1 — s), is 

1 

1 - (m + l)t + st2 -ct^ + 0{t^) ■ 

For any other digraph with s edges and m + 1 vertices, we have that 
the generating series is 

1 

1 - (m + l)t + st2 - dt^ + 0{t^) ' 

with d < c. A perturbation analysis yields that pole of smallest 
modulus is located at + sm~^ — cm~^ in the first case and at 
m~^ + sm~^ — dm~^ in the second, so the first series has smaller radius 
of convergence, hence faster growth of the coefficients. Consequently, 
the first graph has the larger spectral radius. 
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2. The proof 

For any digraph A, let 

oo 

HA{t) = J2xn{A)e, (6) 

n=0 

where Xn{A) denotes the number of walks of length n in A, so that 
Xo(^) = 1; Xi(^) = the number of vertices in A, and Xi(^) = the 
number of edges in A. Let 

R{A) = ^ , (7) 

lim„^oo VXn{A) 

be the radius of convergence of A, and let p{A) = 1/R{A). If the 
adjacency matrix M{A) of A is irreducible then p{A) > is the largest 
eigenvalue of M{A). 

For (m + 1)^ — s > + 1, let VXim^s) denote the finite set of 
digraphs on {l,...,m + l} having precisely (m + 1)^ — s edges. Let 
VVTim, s) C VX{m, s) denote the subset consisting of those digraphs 
whose 0-1 adjacency can be regarded as the Young diagram of a numer- 
ical partition of (m + 1)^ — s; in other words, the rows and columns of 
the adjacency matrix should be weakly decreasing. Then VT>X{m, s) 
is finite, and the cardinality does not depend on m as long as m is 
sufficiently large. Furthermore every digraph in WI{m, s) is con- 
nected in the directed sense, i.e. there is a directed walk between any 
two vertices, digraphs. Hence the adjacency matrix of an element in 
WX{m, s) is by definition irreducible. Furthermore, by a result of 
Schwarz f Sch64j . 

max { p{A) \AeVI} = max { p{A) \ A e VVT{m, s) } . (8) 

Let A denote the complementary graph of A, i.e. the digraph on 
{1, . . . , m + 1} which has an edge i — j iff there isn't an edge i ^ j 
in A. Then the following relation hold (see |iBHV94j and also |Ges77j 

and ^mE\y. 

Hj,{t)H^{-t) = 1 (9) 

If A G VT>X{m, s), then A is a digraph on m + 1 vertices with s edges. 
We have that 

H^{t) = 1 + (m + l)t + st^ + ct^ + 0{t^) (10) 
where c is the number of walks in A of length 2. 

2.1. The case s > 6. Suppose first that s > 6. Backelin |Bac95j 
showed that among all digraphs with s edges, the so-called saturated 
stars have the maximal number of walks of length 2. By a saturated 
star with s = 2k — 1 edges is meant the digraph with a edges (l,i) 
and {i, 1) for 1 < i < k; for s = 2k we add the edge (1, k + 1). So the 
saturated stars with 9 and 10 edges looks as in figure El 
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Figure 2. Saturated star digraphs with 9 and 10 edges 

Note that if i? is a saturated star, then the graph on m + 1 vertices, 
which has an edge i — > j iff (m + 2 — i) — >• (m + 2 — j) is not an edge 
in R, is of the form G{m,s). For instance, the digraphs above have 
adjacency matrices 



11111 
1 
1 
1 
1 



and 



1 1 1 1 1 r 

1 

1 

1 

1 





and if we take m — 7 we get the following adjacency matrices for the 
relabeled complementary graphs: 



■1 


1 


1 


1 


1 


1 


1 


Y 




■1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 
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1 


1 
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1 
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1 
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1 
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1 


1 







1 


1 


1 


1 


1 


1 


1 





1 


1 


1 



















1 


1 





















(12) 



Now suppose that i? is a saturated star, and that A is the digraph 
on m + 1 vertices obtained from R as above. So A and R differ only in 
that A has some isolated vertices. Let i?' be a a different digraph with 
s edges. Let B be the digraph obtained by adjoining isolated vertices, 
so that the total number of vertices becomes m + 1. Let B — B. Then 



Hs{t) = 1 + (m + l)t + st^ + dt^ + o{t^) (13) 

where d is the number of walks in R' of length 2. By Backelin's result, 
c> d. Hence we have that 



1 



1 - (m + l)t + sf^ - ct^ + 0(t4) 

1 



l-(m + l)t + st"^ - dt^ + 0{t^) 
(c - d)t^ + Oit") 



(14) 
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SO we have at once that A has strictly more walks of length 2 than B 
has. By induction, we can show that 

Lemma 1. The exponent off in (fT^ is a polynomial in m of degree 
i — 3, with leading coefficient (i — 2){c — d) . 

Thus, for any j, by taking m sufficiently large, we can achieve 
that the coefficients of f, i < j, in ()14j) are all positive. Recall 
that V'DX{m,s) is finite. Hence, for any j, if we take m sufficiently 
large, then A has the maximal number of walks of length i among the 
B e VVI{m ,s). 

In )Sta97j it is shown that if S5 is a digraph with adjacency matrix 
M, then H^it) = P{t)/Q{t), where Q(t) = det(/ - tM), and P{t) is a 
polynomial of smaller degree then Q. Hence H^(t), H^{t), HA{t), Hsit) 
are all rational functions. Let r = r(m, c) be the pole of HA{t) that is 
closest to origin. Then 1/r is the eigenvalue of M of largest modulus, 
so from the Perron-Frobenius theorem it follows that if is connected 
in the directed sense then r is a positive real number. Writing 

° i-(m + i)i + .t'-ct3 + t^ ,^'';:-';;;s-"- ~ 

\ ' J ' l+alt+a2t^^ hajv2i 2 

we have that r(m, c) is the smallest real root of 

(1 + ait + aat^ + ■ ■ ■ + aN^t^^)^! - (m + l)t + st^ - ct^) 

+ t\bo + bit + --- + bN,t^') = 0. (16) 

We are interested in the asymptotic behavior of r as m — >■ oo. As we 
will show below, we can expand r(m, c) as a Laurent series in m as 

r(m,c) = m~^+sm~^—cm~'^+{2s'^+bo)m~^+(bi—aibo—5cs)m~^+0{m~'^) 

(17) 

Making a similar analysis for the pole r*{m,d) of HB{t) that has the 
smallest modulus, we get 

r(m, c) — r*(m, d) = —{c — d)m^^ + 0{m~^), 

so 

r{m, c) < r*{m, d). 

Thus, for large m, HA{t) has strictly smaller radius of convergence than 
HB{t). When we combine this with Lemma Q we see that by taking 
m sufficiently large, we can achieve that HA{t) ^ Hsit), i.e. that all 
coefficients of HA{t) are > than the corresponding coefficients of Hsit). 
In fact, the inequality is strict for exponents > 2. 



2.1.1. Perturbation analysis of the positive root. We now show how to 
derive the expansion (fTTj) . Replacing m + 1 by 1 /e: in (fTBj) , and clearing 



MAXIMAL SPECTRAL RADIUS OF A DIGRAPH 7 

denominators, we get 

(1 + ait + ast^ H h aN2t^''){e - t + est^ - ect^) 

+ et\bo + bit + ■ ■ ■ + bN,t^') = (18) 

The unperturbed equation is 

(1 + ait + a2t^ + ■■■ + aN^t^^){-t) = 0, (19) 
which has a root at t = 0. We introduce the scahng t = eT and get 

(1 + aieT + a2e^T^ + ■■■ + aN,e^'T^'){e - eT + ehT^ - e^cT^) 

+ e^T\bo + bisT + ■■■ + b^^e'^'T''^) = (20) 

hence 

(1 + aisT + aae^T^ + ■ ■ ■ + aN^e^^T^''){l - T + e^sT^ - e^cT'^) 

+ e^T\bo + bisT + ■■■ + fojVi^'^'T^i) = (21) 

The unperturbed equation is now 1 — T = 0. Hence, T = 0(1) and 
= 0(1), so this is the correct scahng. We make the substitution 
Y = T — 1 and get 

[l+aiEiY+l)+a2e\Y+iy+- ■ .+aj^^e'''iY+lf'] [-Y+ehiY+lf- 

e'ciY + 1)'] + e\Y + l)^(6o + bisiY + 1) + ■ ■ • + bj,,e''' (F + 1)^^) = 

(22) 

It is now clear that Y can be expanded in powers of e, so we make the 
Ansatz 

oo 

Y = J2wie\ (23) 

i=l 

Collecting the coefficients of the powers of e in ()22p we get 

(24) 

-wi (25) 

s — aiWi — W2 (26) 

2swi + ais — c — aiwf — aiW2 — a2Wi — w — 3 (27) 

2sw2 + swf — 3cwi — aic + a2S (28) 

+3aiSWi + &o ~ 2aiWiW2 — aiW3 — 2a2wl — 02^2 — — vl^9) 
These should be zero, which allows us to solve for the Wj's, obtaining 

Wi = 0, W2 = s, W3 = — c, = 2s^ + 60 (30) 

So 

Y = se^ -ce^ + {2s^ + bo)e* + 0{e^), (31) 

hence 

T = 1 + se^ - ce^ + (2^^ + bo)e^ + O(e^), (32) 



e 
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(33) 



This concludes the proof for the case s > 6. 



2.2. The exceptional cases. It remains to take care of the case s < 6. 
Backehn's classification says that if s e {1,3,5} then the saturated 
stars are optimal. Hence, it remains to check s — 2, s = 4, s = 6. 



2.2.1. s — 2. For s — 2 there are two non-isomorphic graphs R, namely 




corresponding to matrices Mi = [H] and M2 = [ ? ] • Recall that 
when Ml or M2 is subtracted from the bottom right corner of the 
(m + 1) X (m + 1) matrix of all ones, the result should be weakly 
decreasing in rows and columns. We call the resulting matrices 



Ai 



111' 



111 
1 



A, 



1 1 r 

1 1 
1 1 



(34) 



The matrix A2 is not weakly decreasing in the last row, so it is not really 
necessary to continue with the calculations, but we proceed anyway in 
order to demonstrate how this is done. We have that 



1 + t 
1-f 



so that 



i^-m+1 



1 — tm — (m — l)i^ 
The smallest positive root of the denominator is 

—m + \/ w? + 4 m — 4 



(35) 



(36) 



so the spectral radius is 



2(m-l) 
2(m- 1) 



—m + V m? + 4 m — 4 



(37) 



(38) 
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2.2.2. s 

digraph 



4. From Backelin's classification we have that for s = 4, the 




with adjacency mat rix [ } } ] has the most walks of length 2. This is 
in accordance with |Fri85j where it is shown that the matrix (jH)) has 
maximal spectral radius among 0-1 matrices with (m + 1)^ — 4 ones. 

2.2.3. s = 6. The remaining exceptional case in Backelin's classifica- 
tion is for s = 6. Then the digraphs in figure El with adjacency matrices 




Figure 3 . Digraphs with 6 edges having maximal num- 
ber of walks of length 2 



Ml 



Mo 



(39) 



1111 
10 
10 


both have 14 of walks of length 2. Note that the first digraph is a 
saturated star. The relabeled complemented matrices are 



1 1 1 
1 1 
1 



1 



11111" 



1 ■■■ 11111 

1 ■•■ 11110 

1 ■■■ 11110 

1 ■■■ 10 



1 



1111" 



1111 
1110 
110 
10 



(40) 



The generating series for the complemented relabeled digraphs are 



l+t-2t^ 



l-mt-mt^ +3t2+2mtf -6*3 



l+2t-t^-t^ 

l+t-mt+3i2-2mt2-2t3+mt3-2t4+mt4 



" P2{t) 



(41) 
(42) 



Regarding the positive root of P2(t) as a function of m, and expanding 
that function as a power series round infinity, we get 



ri(m) -.= 171^ 



m 



7m 



33m-^ + 191m 



-5 



0{m 
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whereas the expansion of the positive root of q2{t) is 

r2(m) := m~^^ — + 7m~^ — 33m~^ + 196m~^ + 0(m~^). 

The first root is therefore shghtly smaller for large m; the difference 
is miniscule, but vive la difference] In fact, since ri(m) — r2(m) < 
for large m, and ri(4) — r2(4) —0.003, it will suffice to show that 
P2it) = 'i2if) = has no solution to demonstrate that ri(m)— r2(m) < 
for m > 4, as shown in Figure |3 Using Macaulay 2 |GSj we can verify 
that {1} is a Grobner bases for the ideal generated by {p2if)i 92 (^)} in 
Q(m)[t]. Hence P2{t),(l2{t) are co-prime in Q(m)[t], so they can not 
have a common zero. 

Figure 4. Difference ri(m) — r2{m) for 4 < m < 10. 








-0.0002 - 
-0.0004 - 




-0.0006 - 




-0.0008 - 




-0.001 - 




-0.0012 - 




-0.0014 - 




-O.0O16 - 




-0.0018 - 




-0.002 - 




-0.0022 - 




-0.0024 - 




-0.0026 - 




-0.0028 - 
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